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Robust Transport Properties in Graphene 
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Two-dimensional Dirac fermions are used to discuss quasiparticles in graphene in the presence of 
impurity scattering. Transport properties are completely dominated by diffusion. This may explain 
why recent experiments did not find weak localization in graphene. The diffusion coefficient of the 
quasiparticles decreases strongly with increasing strength of disorder. Using the Kubo formalism, 
however, we find a robust minimal conductivity that is independent of disorder. This is a conse- 
quence of the fact that the change of the diffusion coefficient is fully compensated by a change of 
the number of delocalized quasiparticle states. 

PACS numbers: 81.05.Uw, 71.55.Ak, 72.10.Bg, 73.20. Jc 
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Recent experimental studies of single graphite layers 
(graphene) have revealed interesting transport proper- 
ties |l|, |2j, I3J]. A quantum Hall effect was found in the 
presence of an external magnetic field with Hall plateaus 
a xy = ±(2n + l)2e 2 /h (n = 0, 1, ...). This result can be 
explained by the band structure of graphene which has 
two nodes (or valleys) due to the hexagonal lattice and 
a linear dispersion around each of these nodes |4| ■ If the 
Fermi energy is near these nodes the quasiparticles can 
be described as fourfold degenerated (two valleys and two 
spins orientations) 2D Dirac fermions. Another interest- 
ing observation is the existence of a minimal conductivity 
cr™ n which occurs if the Fermi energy is exactly at the 
nodes. This quantity shows a remarkable stability with 
a xx n = 3--5 e 2 //i, even if the mobility of the studied sam- 
ples changes by almost a factor 6 from \x = 0.15m 2 /Vs 
to /1 = 0.85m 2 /Vs [l|. The varying mobilities in dif- 
ferent samples indicate a varying density of impurities. 
The transport mechanism must be related to diffusion of 
electrons and holes, caused by impurity scattering. The 
latter can formally be described by a random potential. 
Then, in terms of the Dirac fermions, the impurities can 
create gap fluctuations, i.e. a random Dirac mass, be- 
cause a gap can easily be opened at the nodes of the 
band structure, for instance, by a staggered potential |5(. 

Existing theories of 2D Dirac fermions predicted the 
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value <j xx — e /hit for a single Dirac node 
(i.e. a™" = Ae 2 /hit) such that there is a quantita- 
tive discrepancy by a factor \/it in comparison with 
the experimental observation, as discussed in Ref. |l(. 
fn recent papers several authors applied the Landauer 
formula, instead of the Kubo formula used in previous 
studies, to determine the minimal conductivity also as 
0-™° = 4e 2 [/hit for a rectangular system with aspect ratio 
W/L S> 1 111 H2J . Nomura and MacDonald argued that 
cr™ n could be enhanced by Coulomb scattering, leading 
4e 2 /hit |l3j. Several possibilities for the value 
, using different approaches to the linear response 
were also discussed in Ref. [14]. 

The effect of quantum interference due to impurity 
scattering was studied in a recent experiment [3j. It 
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was found that there is no weak localization in a single 
graphene sheet. Multilayer graphite films, on the other 
hand, exhibit clearly weak localization. These observa- 
tions indicate that graphene has special transport prop- 
erties. 

It will be shown in the following that (1) a calculation, 
based on linear response theory (Kubo formula), gives a 
conductivity of tr™ n = ite 2 /h for the pure system, (2) 
weak scattering leads to a linear Boltzmann conductiv- 
ity similar to what was observed experimentally, and (3) 
there are no weak (anti-)localization corrections due to a 
spontaneously broken supersymmetry which creates dif- 
fusive fermions. 

Starting from the Kubo formula [la ], the conductivity 
tensor a^ v of a system with Hamiltonian H at inverse 
temperature = l/fc^T and for frequency uj reads 
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The current operator is given by the Hamiltonian as 
the commutator j v = —ie[H,r v \. For non-interacting 
fermions with single-particle energy eigenstates \k) (i.e. 
H\k) = ek\k)) a lengthy but straightforward calculation 
yields 

e 2 v-^ 

cv = ^J- z2 ( fc i|[- ff > r Vp2)(&2M£;i) 



ki,k 2 
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with Fermi function fp(e) = 1/(1 + e ,3e ). The identity 
(k 2 \[H,r^]\k 1 ) = (e k . 2 -e kl )(k 2 \r v \k 1 ) and the Dirac delta 
function 7T(5(efc — e) = lim^^o v/(( e k — e) 2 + ?y 2 ) allow us to 
express the conductivity in Eq. (|2|) as a double integral 
with respect to two energies e, e'\ 
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J j ' Tr{[H,r,]6(H - e')[H,r v ]S(H - e)} 
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An alternative expression is obtained by writing the 
Dirac delta functions in terms of the Green's function 
G{z) = (H — z)~ x . Using the identity 

G(z)[H,r^]G(z') = r M G(z')-G(z)r M +G(z)r M (z-z')G(z'), 
the diagonal elements of the conductivity tensor read 
e 2 1 
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Si,S 2 = ±l 
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Me') -Me 



G r 'r{e + is 2 ri2)\- 
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dede',(4) 



where Tr n is the trace related to n additional degrees 
of freedom (e.g., n = 2 for the spinor index of 2D Dirac 
fermions). 

In the special case of graphene the Hamiltonian reads 
in sublattice representation [J, [5[ 



H = a\h\ + a 2 h 2 , 



(5) 



where (Xj (J = 1,2,3) are Pauli matrices. In Fourier rep- 
resentation with wavevector k — (ki,k 2 ) the coefficients 
of the Pauli matrices in a pure system are 

3 3 

hi = — ty cos(Sj ■ k), h 2 = —ty sin(aj • k) 
3=1 3=1 

with the lattice vectors of the honeycomb lattice 2i = 
(-■s/3/2, 1/2), a 2 = (0,-1) and d 3 = (^3/2,1/2). 
iJ can be diagonalized as H = diag(e kl —&k) with 
efe = yhi+h^. The current operator transforms under 
Fourier transformation as — ie[H, r M ] — > edH/dk^. 

There are six points at fc = (±4tt/3\/3, 0), 
(27t/3a/3, ±27r/3), and (— 27r/3v3, ±27r/3) where e^ van- 
ishes, corresponding with the two nodes. The commuta- 
tors are in the diagonal representation of H 
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whose value is 9/4 at all nodes. The 2D k integration at 
each node can be expressed by an integration with respect 
to h\ and h 2 as d 2 k = Jdh\dh 2 , where the Jacobian 
is J = 4/9 at all nodes. Therefore, after the angular 
integration around the nodes the integration is given by 

J(2n/3)e k de k = (8tt /27)e k de k (0 < e k < A). 

This can be inserted in the conductivity of Eq. ([3]) and 
after the summation over all nodes the conductivity 0^ 
reads at low temperatures (/? ~ 00) 
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Inserting the commutators from Eq. ([6]) yields eventually 
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Tr5(2e k - uj)de k 



~k e 
27T 



(7) 



Another factor of 2 comes from the spin-1/2 degeneracy 
of the quasiparticles. Thus our calculation gives for the 
minimal conductivity <7™ n = ire 2 /h. 

In a pure graphene sheet there is only ballistic trans- 
port. Consequently, the diffusion coefficient D is infinite. 
On the other hand, if the Fermi energy is exactly at the 
nodes, the related density of states p vanishes. From this 
point of view, the conductivity, expressed by the Ein- 
stein relation as <J™ n oc pD, depends very sensitively on 
the limits of the model parameters (e.g., the DC limit 
oj — ► 0). A more instructive situation is a system with 
randomly distributed scatterers that may lead to diffu- 
sion (i.e. D < 00) or even to Anderson localization (i.e. 
D = 0) 0, Q3, ES El- A source of disorder in the 
tight-binding Hamiltonian H of Eq. |(SJ) is a randomly 
fluctuating nearest-neighbor hopping rate. For a quali- 
tative discussion of random scattering, the Hamiltonian 
is approximated by the 2D Dirac Hamiltonian 

H D = icriVi + i<J 2 \7 2 + ma 3 . 

A randomly fluctuating gap is introduced by a random 
Dirac mass m with Gaussian distribution of zero mean 
and variance g. The transformation property of the tight- 
binding Hamiltonian 



H 1 = -a 2 Ha 2 



(8) 



is also obeyed by the random Dirac Hamiltonian Hjj af- 
ter rotating ai — > a 2 and o 2 — ► <j\. This property is 
crucial for the formation of a diffusion mode in two di- 
mensions [10|. Models which violate this property, e.g. 
by an additional term proportional to a 2 x 2 unit ma- 
trix, may lead to localization of states near the Fermi 
energy 



18, 1 



Intervalley scattering is ignored by the 
approximation H rs Hd such that we have only inde- 
pendent Dirac cones. The effect of the random mass can 
be studied by applying a perturbation theory, using a 
partial summation of an infinite series of most relevant 
contributions. On the level of the averaged single-particle 
Green's function this leads to a selfenergy term 77 in the 
Green's function: G±(ie) = G(ie ± irj). This is formally 
equivalent to a mean-field approximation of a supersym- 
metric functional-integral approach, where the random 
Dirac mass is replaced by a random supermatrix [lOj. 
Then 77 is obtained as a solution of 



V = 9(V + £ ) / 1(V + e) + fcT kdk/n 



(9) 



The Green's function can be expressed again by a pertur- 
bation expansion, now in terms of fluctuations around the 
mean- field approximation G±(ie). This expansion can be 
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FIG. 1: Conductivity of graphene calculated in mean-field 
approximation (from Eq. (1101) ). 



inserted into the conductivity ([3]). The leading term for 
ie = Ep (Fermi energy) is the Boltzmann conductivity 
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Here a is an approximation of the two commutators in 
Eq. (|3|) and p is the density of states of pure Dirac 
fermions: p(E) = po\E\, where po depends on the cutoff A 
of the spectrum of H. This is a classical result for the con- 
ductivity of Dirac fermions that was already anticipated 
by Fradkin [y] and discussed in the context of d-wave su- 
perconductors by Lee [7| and for graphene by Peres et 
al. 0. The conductivity at E F = is ct™ 11 w 2e 2 ap /h 
and does not depend on r). Thus er™ n is independent of 
the strength of impurity scattering g. Sufficiently away 
from Ep = the conductivity becomes linear, as shown 
in Fig. [TJ This behavior agrees well with the experimen- 
tally observed linear conductivity [l|. 

The next question is whether or not quantum inter- 
ference effects are important. The corresponding cor- 
rections in the conductivity are given by the next order 
terms of the expanded Green's function. This includes 
a logarithmic term (Cooperons) due to a massless mode 
of the fluctuations around the mean-field approximation. 
Previous studies found that the corrections give an an- 
tilocalization effect for conventional scatterers, i.e. an 
increase of the conductivity due to quantum interference 
effects [20] . Additional terms in the Hamiltonian of sec- 
ond order in the momentum can suppress weak antilo- 
calization [2l|. 

A crucial question, implied by the antilocalization ef- 
fect and the absence of weak-localized corrections in ex- 
periments [3(, is whether or not the states in graphene 
are localized. The conductivity in Eq. ((3|) does not di- 
rectly address localization, in contrast to the alternative 
expression in Eq. ((3]) : The long-distance behavior of the 



Green's functions is directly related to the behavior of 
the quantum states. For this purpose we return to Eq. 
(|4]), consider the minimal conductivity (i.e. Ep = 0), 
and take the limits a — > and (3 — ► oo (rji , T)2 — ► are 
implicit) to obtain 
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Assuming that the integrand is finite for < u> <C 1 it 
can be pulled out of the integral for e = 0. The major 
contribution comes from s\ ^ S2, where the poles of the 
Green's function are on different sides of the real axis. A 
finite factor a' takes care of a correction in comparison 
with the exact value of the integral: 

<T « jf^ 2 5>» - r;) 2 rr 2 [G rr .,(c/2 + i Vl ) 



GM-wfi-im)]. (12) 



For localized states the sum, averaged with respect to 
randomness, is finite due to the exponential decay for 
\r — r'\ 3> 1. In the DC limit lo — ► this would lead to a 
vanishing a™ n - 

In order to evaluate the expression in Eq. (fT2f it is 
convenient to return to the functional integral of Rcf. 
[lfj|. It was found that the underlying supersymmetry 
of the integral is spontaneously broken by the mean-field 
solution of Eq. © . The main consequence of this effect 
is a diffusive fermionic mode, similar to the Goldstone 
mode in systems with a rotational symmetry, that can 
be formally described by a complex Grassmann field , F r . 
This allows us to write 



(Tr2[G rr '(ie)G r 'r(-ie)]) = 
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where the action depends on the solution 77 of Eq. © : 
4?7 2 
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g(r) + e) J 4?r(r/ 

The latter contains the diffusion coefficient 
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k 2 ]y k V- k d 2 k. (14) 
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4tt(?7 + e) 



(15) 



that depends strongly on the variance g of the distribu- 
tion of random scatteres. 

The minimal conductivity is obtained from Eq. JT 
for small e = ioj/2, together with Eq. (fj"2"l) , as 
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(16) 



This agrees with the result of the mean-field approxima- 
tion in Eq. (|10p. except for the (undetermined) pref- 
actors. The variation of a' with g for < g < 1 can 



be neglected. Comparing it with the result in Eq. ([7]) 
we conclude that the renormalization factor is a' = n 2 . 
Thus, in contrast to the diffusion coefficient D the mini- 
mal conductivity cr™ n does not depend on g. This implies 
the absence of corrections due to quantum interference 
on characteristic length scales, in agreement with recent 
experimental observations [3J]. 

The fact that <7™ a is so robust with respect to impurity 
scattering can be understood in terms of the Einstein 
relation cr™ n oc pD, where the conductivity is separated 
into the diffusion coefficient D and the averaged density 
of states p at the Fermi energy Ep = 0. The latter is 
calculated from a functional integral similar to Eq. (fl3|) 
as p = rj/ng [lflj. Then with D of Eq. (fT5|) the minimal 
conductivity reads 



In conclusion, we have studied the conductivity of 
graphene, using a model of 2D Dirac fermions. In the 
case of a pure system the ballistic transport leads to a 
minimal conductivity <7™ n = ire 2 /h. In the case of impu- 
rity scattering we found pure diffusion for any strength of 
Gaussian distributed scatterers, where the diffusion coef- 
ficient depends strongly on the distribution. The minimal 
conductivity o-™ n , on the other hand, does not depend 
on the strength of impurity scattering because the change 
of the diffusion coefficient is completely compensated by 
a change of the density of diffusive states. 

This research was supported in part by the Sonder- 
forschungsbereich 484. 
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rj + e r\ + iu/2 



For small e the mean-field equation © gives rj w e n ^ 9 
which implies for D and p 

D « ge^^/An, p « e~ v/g /wg. (17) 

Thus, moving away from the ballistic limit g = 0, the 
conductivity should fall rapidly with increasing random 
potential fluctuations due to a decreasing D in the Ein- 
stein relation. On the other hand, the density of states p 
increases correspondingly so that in er™ n the influence of 
random scattering is compensated. This is in agreement 
with the direct evaluation of ct™ 11 in Eq. (fl6)l . The re- 
sults for p and D in Eq. (|17[) describe a nonperturbative 
effect of disorder which is not visible within an expansion 
in powers of g. 

Strong potential scattering by charged impurities in 
the substrate, for instance, can lead to a destruction 
of the massless fermion mode used in Eq. (fi"4|) . This 
can cause localization and a vanishing cr^ 1 ", at least at 
very low temperatures. The localized regime cannot be 
treated within a conventional field theory but would re- 
quire either a numerical finite-size scaling approach [22j 
or a strong-disorder expansion. 
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